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Abstract 



We assume the validity of the Standard Model up to an arbitrary high-energy scale 
and discuss what information on the early stages of the Universe can be extracted from 
a measurement of the Higgs mass. For Mh ~ 130 GeV, the Higgs potential can develop 
an instability at large field values. From the absence of excessive thermal Higgs field 
fluctuations we derive a bound on the reheat temperature after inflation as a function 
of the Higgs and top masses. Then we discuss the interplay between the quantum 
Higgs fluctuations generated during the primordial stage of inflation and the cosmo- 
logical perturbations, in the context of landscape scenarios in which the inflationary 
parameters scan. We show that, within the large-field models of inflation, it is highly 
improbable to obtain the observed cosmological perturbations in a Universe with a light 
Higgs. Moreover, independently of the inflationary model, the detection of primordial 
tensor perturbations through the 5-mode of CMB polarization and the discovery of a 
light Higgs can simultaneously occur only with exponentially small probability, unless 
there is new physics beyond the Standard Model. 



1 Introduction 



The search for the Higgs boson and the measurement of its properties are one of the primary 
goals of the LHC. The mere discovery of a Higgs can be viewed as a possible indication of 
additional new physics not far from the electroweak scale, because of the high sensitivity 
to short-distance quantum corrections of the mass term associated to a fundamental scalar. 
However, even in the absence of any new-physics discovery at the LHC, a measurement of the 
Higgs mass can give us useful hints on the structure of the theory at very high energies. This 
is because, at large field values, the Higgs potential can develop an instability or become non- 
perturbative, depending on the precise value of the Higgs quartic coupling A or, ultimately, 
on the Higgs mass Mh- Because of the logarithmic dependence of A on the energy, such 
considerations [H El El HI E] can test the properties of the theory up to extremely small 
distances, which are otherwise totally unaccessible to any imaginable collider experiment. 

In this paper, we want to use the same considerations for a different purpose. Rather 
than trying to infer new particle-physics properties at small distances, we will assume the 
validity of the Standard Model (SM) up to an arbitrary high-energy scale, and find what 
information on the early stages of the Universe can be extracted from a measurement of the 
Higgs mass. We will first obtain that in the Higgs mass range 114 GeV < < 130 GeV, 
where the electroweak vacuum is potentially metastable, the absence of excessive thermal 
Higgs field fluctuations in the early Universe imposes a bound on the reheat temperature 
after inflation Tr#. 

Then we will discuss the interplay between the quantum Higgs fluctuations generated 
during the primordial stage of inflation and the cosmological perturbations which are either 
currently observed in the form of Cosmic Microwave Background (CMB) anisotropics, or 
might be detected in the near future in the form of tensor (gravity waves) perturbations. The 
key ingredient is that all these fluctuations depend upon the Hubble rate during inflation and 
therefore, under certain assumptions, it is possible to relate the amount of Higgs fluctuations 
to observable properties of the CMB. However, excessive fluctuations of the Higgs field during 
inflation can pose a threat to the stability of the present vacuum, if the Higgs mass lies in 
the metastability window. 

We will assume that the various initial inflationary patches of the Universe are character- 
ized by different microphysical parameters. In this sense we take the point of view that the 
underlying theory has many vacua, realized in different patches of the Universe. This pic- 
ture, usually referred to as the landscape [6], has been put forward especially in the context 
of string theory. Under this assumption, from CMB observations and from a measurement 
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of the Higgs mass, we can derive probabilistic conclusions on the properties of our Universe. 

In particular, within the class of large-field models of inflation, we will compute the 
probability to have a Universe at the end of inflation which both survived the quantum 
Higgs fluctuations and has the right amount of observed cosmological perturbations. Such 
probability is extremely (exponentially) small, if the Higgs mass is below 124 GeV (for the 
present central value of the top mass). Moreover, we find that the discovery of a light Higgs 
together with the detection of primordial tensor perturbations through the 5-mode of CMB 
polarization (at a level quantitatively described in sect. 6) would imply that we live in a very 
atypical Universe, whose probability decreases exponentially when the Higgs mass decreases. 
Such discovery could be interpreted as an indirect evidence for the existence of new physics 
beyond the Standard Model at some intermediate energy scale. 

The paper is organized as follows. In sect. 2 we briefly review the Higgs mass instability 
window. In sect. 3 the bounds on the reheating temperature after inflation are discussed. 
In sect. 4 we compute the survival probability of the electroweak vacuum during inflation 
and in sects. 5 and 6 we relate it to the curvature and tensor perturbations, respectively. 
Section 7 states our conclusions. The appendix contains technical details relevant to the 
calculation of the Higgs mass instability window. 



2 The Higgs mass instability window 

Let us start by reviewing the range of Higgs masses for which the electroweak vacuum is 
metastable. Since our considerations refer to large field values, for our purposes it is perfectly 
adequate to neglect the bilinear term and to approximate the potential of the real Higgs field 

h as 

V(h) = (1) 

We will work in next-to-leading order approximation and include the field-dependence of the 
quartic coupling A, determined by the two-loop renormalization-group (RG) equations, as 
well as the one-loop corrections to the effective potential (as in ref. [3]). The couplings A 
and the top-quark Yukawa coupling h t , entering the RG evolution, are then related to the 
physical Higgs and top pole masses. The explicit expressions are given in the appendix. 

The request that the electroweak vacuum (h) = v, with v = {^/2G F )" 1 ' 2 = 246.22 GeV, 
is the true minimum of the potential, up to a cut-off scale A, implies A(/i) > for any fi < A. 
This condition is not satisfied at some energy scale whenever Mh is below some critical value 
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Figure 1: The instability scale A as a function of the Higgs mass for three different values 
of the top mass M t . 

M£ given by 

* < * - ^« (^^)— V (^^) ±2 GeV. 

(2) 

We have explicitly shown the dependence on the two most important SM parameters, nor- 
malizing their effects in units of one standard deviation from their experimental central value, 
taking M t = 170.9 GeV ± 1.8 GeV [7J and a s (M z ) = 0.1176 ± 0.0020 [g]. Besides the uncer- 
tainties in the SM parameters, eq. ([2]) has an overall error due to higher-order corrections 
which, parametrically, are expected to shift the result by an amount 0(a 2 s Mij7t 2 ). However, 
since the two-loop QCD correction to the top-quark pole mass (included in our calculation 
and shown in the appendix) has a large coefficient, we conservatively estimate the theoretical 
error to be 2 GeV. Figure [T] shows the instability scale A, at which the quartic coupling A 
becomes negative, as a function of the Higgs mass for three different values of the top mass. 
At scales larger than A, the Higgs potential becomes negative, and then it develops a new 
minimum. The endpoints of the lines, marked with a square, correspond to M£ and to the 
scale at which the new minimum of the Higgs potential characterized by a large vacuum 
expectation value becomes degenerate with the electroweak one. 

A lower bound on Mh is obtained by considering the instability of the electroweak vacuum 
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Figure 2: Lower bounds on from absolute stability (upper curves) and T = metasta- 
bility (lower curves). The width corresponds to a s [Mz) = 0.1176 ± 0.0020 (with the higher 
curve corresponding to lower a s ) and we do not show the uncertainty from higher-order 
effects, which we estimate to be below 2-3 GeV. The horizontal line is the LEP mass bound. 
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under quantum tunneling. Since the vacuum transition is dominated by late times, this 
bound is independent of the early history of the Universe. The tunneling probability p is 
given by 

p = max Vuh^exp ( — ) . (3) 



h<A ru,u ^ v 3|A(/i)i; • 

We have also included one-loop corrections to the bounce configuration [2], which however 
have only a small impact on the final bound on Mh, contributing to less than a GeV. In 
eq. (j3j), Vu is the space-time volume of the past light cone of the observable Universe, and 
we take Vu = t v , where tu is the lifetime of the Universe. Taking tu = 13.7 ± 0.2 Gyrs from 
WMAP data [9], the metastability limit p < 1 imposes the bound on the Higgs mass [2] 

M h > 105.6 GeV + 5.2 GeV f M '- 17 " 9GeV V2.5 GeV ( = °' 1176 WeV. (4) 

V 1.8 GeV J V 0.0020 J V ; 

The error of 3 GeV is estimated by combining uncertainties from higher-order corrections 
and from the prefactor in p. 

The bound in eq. (HI) is useful only when it is stronger than the direct experimental limit 
on the Higgs mass [TO] 

M h > 114.4 GeV at 95% CL. (5) 

In summary, eqs. (jSj), (j3J) and (jSJ) define the Higgs-mass window (see fig. [2]) in which the elec- 
troweak vacuum is metastable and therefore potentially sensitive to large field fluctuations 
during the early stages of the Universe. 



3 Bound on the reheating temperature 

At sufficiently large temperatures, thermal fluctuations in the early Universe plasma can 
trigger the decay of the metastable electroweak vacuum [12l HI [13] by nucleation of bubbles 
that probe the Higgs instability region. On the other hand, high-temperature effects also 
modify the Higgs effective potential, with a tendency of making the origin more stable. 
The contribution of the different plasma species to the potential (or rather, free energy) 
in the non-interacting gas approximation is given by standard one-loop (bosonic/fermionic) 
thermal integrals. Each particle species, with /i-dependent mass M a (h), contributes to the 
free-energy 

e -y/x*+Mg(h)/T*~ 

+ ^^^{M^-KW + H^T 2 )] 372 } , (6) 



6 a V{h) 



N a e a I dx x 2 log 
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where N a counts the number of degrees of freedom, e a = +1(— 1) for bosons (fermions) and 
U a (h,T 2 ) is the thermally corrected mass of the corresponding species (see e.g. ref. |14j). 
The second line in eq. (|6]) takes into account the effect of resumming hard-thermal loops for 
Matsubara zero modes. For our numerical work we used a series expansion of these integrals 
in terms of modified Bessel functions [15], avoiding high T expansions. 

The energy E C (T) of the smallest critical bubble large enough to grow (overcoming the 
surface tension penalty) controls the false vacuum decay rate through a Boltzmann sup- 
pression factor exp [—E C (T)/T]. The quantity E C (T) is computed by solving for the 0(3) 
bounce solution [12] using the finite T potential described above. It is easy to show [I] that, 
parametrically, E C {T)/T ~ 7rg/|A(T)|. 



The vacuum decay rate per unit volume is 

r(T) ~ t 4 r^^ 1 



exp[-£ c (T)/T] . (7) 



2vrT 

The differential decay probability dP/dlnT is obtained by multiplying T(T) above by the 
volume of the Universe at temperature T and the time spent at that T. In a radiation 
dominated Universe one has 

dlnT- L[I)Tu T* \T 

where T ~ 2.73° K ~ 2.35 x 1(T 4 eV and M p = 1.2 x 10 19 GeV is the Planck mass. The 
previous result assumes T is smaller than the reheating temperature after inflation, Tjuj. 
For temperatures T > Trh, the previous result gets modified to 

dP WrT7 , ,M P (T RH \ W ( T x:i 



d]nT -nT)r^[^) , (9) 

to take into account the period of expansion during inflaton dominance before reheating is 
completed. The final decay probability P, resulting from integrating in dlnT the previous 
expressions, is not a properly normalized probability. Its interpretation is that the fraction 
of space converted to the true vacuum goes like e~ p . 

For a given value of the Higgs mass in the metastability window of fig. [2] the requirement 
that the false vacuum does not decay during the high T stages of the early Universe (that 
is, e~ p < 1) will set an upper bound on Trh- In fact Trh is not the maximal temperature 
achieved after inflation. Such maximal temperature occurs after inflation ends and before 
reheating completes and is given by [TB] 

/ Q \ 2/5 / r \ 1/8 1/8 (rp \ 

J- max V 8 / U 3 J n^tT \ P f ^' ^ } 

\ / n / i/* y-Lmax) 
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Figure 3: Upper bounds on T RH} as functions of M h , from sufficient stability of the elec- 
troweak vacuum against thermal fluctuations in the hot early Universe for three different 
values of the top mass. The lower curves are for Hf = 10 13 GeV, the upper ones for Hf 
deduced from eq. (jTTj) . Hf = [A7r 3 g^(TjiH)/ / i5] 1 ^ 2 (T RH /M p ), which corresponds to the case 
of instant reheating. We take as(M z ) = 0.1176. Lowering (increasing) as(M z ) by one 
standard deviation lowers (increases) the bound on T RH by up to one order of magnitude. 
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where g*{T) counts the effective number of degrees of freedom (with a 7/8 prefactor for 
fermions) with masses <C T and Hf is the Hubble parameter at the end of inflation. The 
metastability bound on Trh therefore depends on the particular value of Hf-. for a given 
Trh, the value of T max grows with Hf. Therefore the metastability constraint on Trh will 
be more stringent for larger values of Hf. 

Figure [3] shows this metastability bound on Trh as a function of the Higgs mass for 
various values of the top mass and for two choices of the Hubble rate Hf at the end of 
inflation. The lower curves correspond to Hf = 10 13 GeV while the upper ones have Hf = 
[4tt 3 #* (T rh ) /45] 1/2 (T% H /M P ), the lowest value allowed once it is required that the inflaton 
energy density = 3M^Hf /(8tt) is larger than the energy density of a thermal bath with 
temperature T RH 

45 



T RH < 



lA Ml^H)' 2 . (11) 



^g*{T RH ) _ 

This condition also corresponds to the requirement that the inflaton lifetime T^ 1 is larger 
than the Hubble time H J . The bound on Trh gets weaker for smaller values of the top mass 
(or larger values of the Higgs mass) since the instability scale becomes higher and eventually, 
for P3] 

A4 > 117.4GeV + 4.2GeV ( M > ~ 17 ^ 9 °» V ) - 1.6G.V ±3 GeV, 

V 1-8 GeV J V 0.0020 J 

(12) 

the bound on Trh is lost: the vacuum is sufficiently long-lived even for T ~ M p [that is, 
J^ p (dP/d In T)d In T < 1]. This is the reason why the lines in fig. [3] stop at some value of 
M/j. From this figure it is also clear that the bound on Trh is very sensitive on the value of 
M t , with the experimental error in M t being the main source of uncertainty. 

It is natural to ask what implications the bound on Trh shown in fig. [3] has for leptogen- 
esis. SM thermal leptogenesis, with hierarchical right-handed neutrinos, sets a lower bound 
on Trh as a function of Mi, the mass of the lightest right-handed neutrino [17] . This bound 
reaches its minimum for M x ~ Trh, when Trh > 3 x 10 9 GeV [18] . This condition could 
be in conflict with the upper bound on Trh shown in fig. [3l if the Higgs mass turns out to 
be very close to the LEP lower limit and if the top mass is on the high side of the allowed 
experimental range. However we stress that these considerations apply only to the case of 
hierarchical thermal leptogenesis in the SM, with no new physics present below the scale 
Mi. 

The Yukawa couplings h u of the heavy right-handed neutrinos could in principle affect 
the bound on T RH , since h v can modify the instability scale of the Higgs potential [19j with 
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its effect on the evolution of A above the M\ threshold. Because h 2 u = m u Mi/v 2 , such effects 
turn out to be important only if the mass of the right-handed neutrinos is sufficiently large, 
Mi > (10 13 — 10 14 ) GeV [19|. Therefore, the existence of heavy right-handed neutrinos 
can modify the bounds on T RH we have obtained only at such large energy scales, i.e. for 
Trh > Mi > (10 13 - 10 14 ) GeV. 

4 Survival probability of the electroweak vacuum dur- 
ing inflation 

In the previous section we have discussed the stability of the electroweak vacuum against 
thermal fluctuations. These are expected to drive the Higgs field towards the instability 
region if Trh is larger than the critical temperature shown in fig. [3j If reheating is an 
instantaneous process, so that Hf ~ g X J 2 (Trh)(T rh / M p ) , then the upper bound on Trh 
(which in this case coincides with T max ) can be translated into an upper bound on Hf. We 
obtain that (for M t = 170.9 GeV) H f /A < 0.1 for M h = 115 GeV and H f /A < 10 4 for 
Mh = 117 GeV. Not only this limit on Hf becomes quickly weaker as Mh is increased inside 
the metastability window, but it also completely evaporates if reheating is a prolonged 
process with T max 3> Trh- Moreover, during inflation the Hubble rate could be (much) 
larger than its value at the end of inflation. Therefore, even if thermal fluctuations do not 
destabilize the electroweak vacuum, the ratio Hj A during inflation can be much larger than 
one. It is thus interesting to address the issue of what is the fate of the electroweak vacuum 
when we account for the Higgs fluctuations generated during inflation, whose amplitude is 
directly proportional to the Hubble rate and can probe the instability region. 

In the inflationary cosmology picture, the present structure in the Universe is supposed 
to originate from vacuum fluctuations which were quantum-mechanically excited during an 
inflationary stage [201 [21] . Indeed, any scalar field whose mass is lighter than H, where H 
denotes the Hubble rate during inflation, gives rise to an almost scale-invariant spectrum 
of perturbations on superhorizon scales [22l [23]. The Higgs field, if light enough during 
inflation, does not represent an exception. We will start by considering the case in which the 
Higgs field is minimally coupled to gravity and therefore it can be considered nearly massless 
during inflation, as long as H ^> Mh- At the end of sect. Owe will discuss the interesting 
case in which there is a direct coupling between the Higgs bilinear and the Ricci scalar. 

If the Higgs mass lies in the metastable window (114 GeV < Mh < 130 GeV), excessive 
fluctuations of the Higgs field may pose a threat to the stability of the electroweak vacuum. 
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In the inflationary picture, long wavelength perturbations of the Higgs field may be gener- 
ated and they behave in the same way as a homogeneous classical field. We suppose that 
in the beginning of inflation the whole inflationary domains are characterized by a very ho- 
mogeneous Higgs field with vanishing vacuum expectation value. Then the domains become 
exponentially large. The Higgs quantum fluctuations divide the Universe into exponentially 
large regions with different values of the Higgs field: the Universe becomes quickly filled with 
domains in which h typically changes from — (h 2 ) 1 ^ 2 to (h 2 ) 1 ^ 2 . Here (h 2 ) is the variance of 
the Higgs field which we will estimate in the following. 

The process of generating a classical Higgs field configuration in the inflationary Universe 
can be interpreted as the result of the Brownian motion of the Higgs scalar field under the 
action of its quantum fluctuations which are converted into the classical field when their 
wavelengths overcome the horizon length. The best way to describe the structure of the 
Higgs fluctuations is provided by the stochastic approach in which one defines the distribution 
of probability P c (h,t) to find the Higgs field value h at a given time at a given point [24|. 
Perhaps a more adequate approach might be based on the distribution of probability P p (h, t) 
to find the Higgs field value h at a given time in a given physical volume [25] which takes into 
account the exponential growth of the volume of the various domains. Both distributions of 
probabilities are plagued by various problems (for a recent discussion on this issue related 
to the idea of eternal inflation, see ref. [26]). The co moving probability P c at a given point 
ignores the possible creation of new volumes (new points) during inflation; the physical 
probability suffers from normalizability problems. In this paper we infer our results by 
making use of the comoving probabiliy P c which is technically more feasible. The subscript 
c serves to indicate that P c corresponds to the fraction of original comoving volume filled 
by the Higgs field h at the time t. The comoving probability satisfies the Fokker-Planck 
equation 

dP c _ d 
~dt ~ dh 

In writing this equation we have supposed that the Higgs field gives a negligible contribution 
to the energy density of inflation and that the probability does not sensitively depend upon 
the value of the inflaton field. The first assumption always holds in the relevant region of 
Higgs field configurations and we will come back to the second point later on. 

To solve eq. f|T3|) we assume that the Hubble rate is approximately constant during infla- 
tion. We denote from now on by A the value of the Higgs field at which the maximum of the 
potential V(h) occurs. It is very close to the value of the energy at which the Higgs potential 
vanishes shown in fig. [TJ In those patches where the Higgs field takes any value larger than 
A, the Higgs field rapidly rolls down towards the region where the potential is negative, 



H 6 dP, 



8tt 2 dh 



3H 



(13) 



10 



V(h) = — \X(h) \ h A /4. Therefore, the domain where the barrier has been surmounted will 
correspond to bubbles of AdS space which will rapidly collapse and disappear. To account 
for this process we impose that the solution of eq. ( fl3l) satisfies the boundary condition 



P c (A,t) = 0. (14) 

Physically, this condition dictates that, once the unstable region of the Higgs potential is 
probed, the probability to jump back to the stable region vanishes. As dP c /dh does not 
vanish at h = A, the probability that the Higgs field is in the stable region h < A will 
decrease with time. Since P c (h, t) gives the fraction of the comoving volume occupied by 
regions inside which the Higgs field takes the value h at a given time t, we expect that this 
fraction will become smaller and smaller as time goes by during the inflationary stage. 

To investigate the amount of Higgs fluctuations during inflation, we can also parametrize 
them in terms of the Higgs correlations determined with the help of the comoving probability 
given by eq. (|T3j) . 

d H 3 m H 3 \ m 

Jt (h m ) = — 2 m( m ~ l )^ hm ' 2 ) - ^^V'(h)) + ——^-P' c (h) , (15) 

where m is an even integer, the prime denotes derivatives with respect to h, and 

(h m ) = [ dh h m P c (h,t) . (16) 

J -A 

We consider the case in which P c (h, t) is an even function of h, so that any Higgs correlation 
functions with odd m identically vanishes. 

For constant H, eq. ( fl3l) can be solved by separation of variables 

P c (h,t) = f^c n e< aV+3 ^ 1 )^ n (h), a = ^ , (17) 

n=0 

where $ n and a n are the eigenfunctions and the eigenvalues of the equation 

< - = -aA . (is) 

Since V(h) is an even function of h and eq. ( Tl8l) does not mix even and odd eigenfunctions, 
we look for solutions which are even, in the range —A < h < A. Moreover, since the condition 
in eq. (fl4j) has to be satisfied at any time t, we have to impose $ n (A) = 0, for any n. 

We assume that the Higgs field is initially localized at h = 0, P c (h, 0) = 5(h), and study 
its evolution. We have solved the Fokker-Planck equation numerically, but it is useful to 
give here approximate analytic solutions in order to describe our results. 
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(i) Case H ^> A. Let us first consider the most interesting case in which the Hubble 
rate is much larger than A. In this case the potential term in eq. (1181) can be neglected since 
aV'$' n /$n < A(A/if) 4 «C 1, once we use the condition \h\ < A. This means that the loss of 
probability caused by the AdS instability dominates the dynamics, while the effects from the 
Higgs potential are negligible. Then the solution of the Fokker-Planck equation satisfying 
the appropriate boundary conditions k *] 



2 ) 8A 7 COS 



n=0 



l\ 7lh' 

" + 2 X 



(19) 



Notice that, in the limit of large A, the probability distribution in eq. (fT9l) reduces to a 
Gaussian 

/ 2tt 2A 2 , . 

as it can be easily obtained by switching from discrete to continuous variables (n+l/2)7r/A — > 
A; and by integrating over fc. 

The survival probability Pa for the Higgs to remain in the region | h\ < A and the variance 
of the Higgs field are given by 



P A (t)= f\hP c (h,t) = -f]^ T e^) 2 ^, 



(h 2 ) 



2A 2 



dh h 2 P c (h,t) = — V 

7T f ^ 



n=0 ^ 



At small t, eqs. (J2TJ) and (J22J) become 

P A (t) ~ 1 - 
</> 2 > 



P 3 t 
2tt 3 A 2 

Hh 



e H 6 t 



An 2 



(n+±)V 

(small t) 
(small t) 



2 I 8A 2 



(21) 



(22) 



(23) 
(24) 



At the initial stages, the loss of probability is very small and Pa is exponentially close to 
one, while (h 2 ) starts growing linearly with time. To study how fast the probability decays 
with time, we can consider eqs. ( 1211) and (1221) in the limit of large t, 



P A (t) ~ — e 32A 2 " 
7T 

2v _ 4(tt 2 -8)A 2 ^ # 



(h 2 ) 



e 32A 2 



7T" 



[H > A, large t) 
[H > A, large t) 



(25) 
(26) 



1 One can show that P c (h, 0) = 5(h) using J2^=o cos[(n + 1/2)^] = tt8(x). 
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The survival probability Pa and the Higgs correlation function exponentially decay with 
time. The normalized Higgs variance, which measures the correlation function within the 
surviving domains, tends asymptotically to a constant, (h 2 )/P\ — > (1 — 8/7r 2 )A 2 ~ (0.44 A) 2 . 
The asymptotic behavior shown in eqs. (I2T!) and (|22|) is valid for t ^> A 2 /H 3 and therefore it 
is reached very rapidly, justifying the assumption of neglecting the time dependence of the 
Hubble constant. 

From eq. (1221) we can also obtain the maximum value of the Higgs variance 

(/i 2 ) max ~ (0.39 A) 2 , (27) 

which is obtained at the time t ~ 10.6 A 2 /H 3 . 

(ii) Case A 3> H. At the first stages of inflation, the first term in eq. (fT8l) dominates and 
the results given in eqs. (1231) and ( 1241) for small t are valid also in the case A ^> H. At later 
times the shape of the probability distribution changes because the second term in eq. ( |T8l) 
becomes relevant. To estimate when the effect of the Higgs potential is important, we study 
eq. ffT5l) for the variance of the Higgs field (m = 2) 



SC» 2 ) = ^[^ + A^(A)]-|(AV. (2f 



Here we have neglected the contribution to the correlators from the Higgs dependence in the 
Higgs quartic coupling, and we have adopted a Hartree-Fock approximation for the nonlinear 
term by taking {hV'(h)) ~ 3X(h 2 ) 2 , where A is evaluated at the scale a/ (h 2 ). From eq. ( 1281) 
we deduce that the Higgs variance starts growing linearly with time, (h 2 ) = H 3 t/(4ir 2 ), and 
then it reaches a maximum value which can be estimated by equating the right hand side of 
eq. ( 1281) to zero 

H 2 

(/i 2 ) max ~ — . (29) 



This value is reached at a number of e- foldings approximately given by N ~ Ht ~ 7Ta/2/A. 
Therefore this maximum is reached very promptly which, also in this case, justifies neglecting 
the time dependence of the Hubble rate. At times t > 1/ (\f~\H), the friction term starts being 
relevant in eq. ( JTBl) . In the absence of the AdS instability region, the comoving probability 
would reach a stationary form P c (h) oc exp[— 87r 2 V(h)/3H 4 }, i.e. ao = and $o — 1 i n 
eq. (|T7|) . However, in the presence of the AdS region, this solution has to be modified to 
account for the loss of the probability. The stationary solution does not satisfy the boundary 
condition in eq. (fT4"|) but, for A ^> H, P C {A) is exponentially close to zero. Therefore we can 
obtain the correct solution by perturbing around the stationary solution valid for A — > oo, 
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and we find $q(/i) ~ 1 - 1(h) /I (A) and a ~ I /1(h), where 

J(/i) = [ k dh'e- aV ^ f dh"e aV ^ . (30) 

JO Jh' 

At late times, the survival probability decays with time as P\(t) oc e -7 *, with 

1 = ^"f(T7i)(T) (A » B, large t) . (31) 

Here we have integrated eq. (1301 neglecting the field dependence of the quartic Higgs coupling 
A and we have taken the limit A 3> H. In the more relevant case in which the Hubble rate 
is changing with time, we have numerically checked that the survival probability decays as 
exp[— f* dt /r y(t')], where one can approximately use the expression of 7 in eq. (I3T]) simply 
plugging the appropriate time-dependence of the Hubble rate. This is also true in the case 
H A, using the time dependence shown in eq. (j2"5|) . 

5 The survival probability and the comoving curvature 
perturbation 

From the results of the previous section we deduce that, if H ^> A, an exponentially small 
fraction of the initial comoving volume survives the Higgs quantum fluctuations during infla- 
tion. The number of surviving domains will be suppressed by exp (—H 2 N/32A 2 ), where N is 
the total number of e-folds. May we conclude that our Universe is very unlikely? The answer 
is no. Suppose, indeed, that to quantify the probability of ending up in a "well-behaving" 
domain (the vacuum has the correct Higgs vacuum expectation value) we simply count the 
fraction of domains where the value of the Higgs field is less than A. This fraction is given 
by the comoving survival probability. However, those regions where the value of the Higgs 
field has become larger than the critical value during the inflationary stage simply do not 
exist. As soon as the instability point is reached, those regions collapse and disappear. Only 
the "well-behaving" regions remain and their corresponding fraction is unity. 

Can we then conclude that the parameter space where Pa is exponentially small is ac- 
ceptable? As we will show, the answer is again no, at least in a probabilistic sense. Indeed, 
in the class of inflationary models dubbed large-field models, the "well-behaving" domains, 
where the vacuum expectation value of the Higgs is correct, typically have an insufficient 
amount of curvature perturbations generated during inflation. In the inflating domains that 
are characterized by the correct amount of curvature perturbations, the Hubble rate, and 
therefore the square root of the Higgs variance, turns out to be much larger than the scale 
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A. Therefore, among the surviving regions, the ones having the right amount of curvature 
perturbations will be exponentially rare. 

Observations of the CMB anisotropies are consistent with a smooth and nearly Gaussian 
power spectrum of curvature perturbations with an amplitude 



Cobs ~ 5 x 10 

Inflation predicts such a spectrum with amplitude 

m 



-5 



c 



2tt( 



(32) 



(33) 



where if* and 0* are the value of the Hubble rate and the inflaton field, respectively, when 
there are about 60 e-folds to go till the end of inflation. The dot stands for differentation 
with respect to time. 



Let us focus on the large-field models of inflation. We parametrize their potential by 

A-p 



v{<j>) 



V 



(34) 



with p a positive integer. This determines the amplitude of the curvature perturbations 



c 



4iV t 



1/2 



P* 



(35) 



irp y Mp 

where N* is the number of e-folds till the end of inflation when the observable scales leave the 
horizon^. Notice, in particular, that the perturbation is directly proportional to the Hubble 
rate. This will play a crucial role in the following because larger values of the Hubble rate, 
and therefore larger values of (, tend to destabilize the electroweak vacuum. 



In the case in which H ^> A, the survival probability becomes 

H?N 



Pa oc exp 



o 



32A 2 



exp 



16(p + 2)A 2 



(36) 



where Hi is the Hubble rate at the beginning of inflation and N is the total number of e-folds. 
Using eq. fl35|) and (Pj/P*) = (iV/iV*)^ 4 , the survival probability can be expressed in terms 
of the total number of e-folds and the comoving curvature perturbation 

2 . ... E±ll 



Pa oc exp 



irp ( (M, 



p + 2 V 8A 



N 



(37) 



2 The number of e-folds to go till the end of inflation for those scales which exit the horizon and are 
today relevant for observation depends on the reheating temperature Trh through the relation N* ~ 60 + 
1/6 ln(— tit) + l/31n(Tfjjj/10 GeV) — 1/3 In 7 [21], where nr is the spectral index of tensor perturbations 
and 7 is the ratio of the entropy per comoving volume today to that after reheating. 
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We now take a crucial step and assume that initially the Universe is characterized by 
local domains where the inflationary parameters may assume any value. As mentioned in 
the introduction, this assumption is inspired by the string landscape picture [6], in which 
the theory has an enormous number of vacua, each determining different values for the un- 
derlying parameters. This hypothesis agrees especially well with the point of view followed 
in this paper of ignoring the hierarchy problem associated to the Higgs mass and extending 
the validity of the SM up to very high-energy scales. Indeed, this approach could be justified 
in the context of a landscape scenario in which the Higgs mass parameter scans among the 
different vacua. In practice, in the case under consideration, we allow for the possibility that 
the various initial comoving patches of the Universe are characterized by different micro- 
physical parameters, such as the parameter \i in eq. (|34p and the initial value of the inflaton 
field 0j. These two parameters determine the other cosmological variables 



P+2 4-p 
2 



Using eq. fl38|) . we trade the parameters fa and /i for the two more physical quantities N and 
(. The way the inflationary parameters are distributed at the beginning of inflation among 
these different regions is unknown and we take it to be simply flat in N and (. Our results 
are not very sensitive to this assumption (as long as the parameter distributions are not 
sharply peaked), because the probability function we are interested in has an exponential 
dependence on our variables. Among the comoving patches which have survived the Higgs 
instability, the probability that the total number of e-folds is between N and N + dN and 
that the comoving curvature perturbation is between ( and ( + d( is Pa{(, N) d( dN, with 
Pa(C) N) given in eq. ([3] 



The range in which the parameters scan is in principle arbitrary. We can restrict it using 
anthropic priors, although this is not essential for our conclusions. If perturbations are too 
large, ( > 10~ 4 , they become non-linear when the average density of the Universe is too 
high to allow the resulting structures to guarantee a stable environment for life |27J. On the 
opposite side, if the perturbations are too small, ( < 10~ 6 , the majority of the overdense 
regions are not capable of cooling quickly enough to form structures (27] . Therefore, we scan 
the comoving curvature perturbations only over the anthropically allowed region 10~ 6 < ( < 
10" 4 . 



The probability in eq. ( 137|) is dominated by the smallest value of N, which we choose to 
be iV* in order to allow at least sixty e-folds of inflation. Integrating over the total number 
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Figure 4: The quantity log 
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as a function of the Higgs mass 



of e-folds, the probability becomes 



P A (0 <x E^_(a) 



a: 



Tip f(M } 



(39) 



a p + 2 V 8A 

Here E n (x) is the exponential integral function and we have expanded for large values of 
the parameter a, which is appropriate since eq. ( 1371) is valid only for H ^> A (i.e. a ^> 1). 
When a < 1, -Pa(C) should be replaced by a constant. The probability -Pa(C) is peaked at 
the smallest value of the cosmological perturbation (. Therefore, the fraction of domains 
whose value of the cosmological perturbation is as large as the measured Cobs is extremely 
small, even if we impose the prior of restricting our considerations only to values of ( within 
the narrow anthropic range: 



10 



-6\ 3 



exp 



8A 



(40) 



J c>10 -6^C^a(C) V Cobs J v 

We recall that eq. (T4"0l is valid for H ^> A which holds in the anthropic window for ( as 
long as A < 10 12 GeV. In fig. ([5]) we plot this ratio as a function of the Higgs mass and for 
different values of the top mass and for p = 2. As expected, if the Hubble rate is large 
enough to account for (, see eq. (1351) . the fraction of volumes which survived the quantum 
Higgs fluctuations is extremely small, unless the mass of the Higgs is close to (or larger than) 
the critical of eq. ([2]) shown by the upper curves in fig. 2. 
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Let us now discuss how our results change if we adopt a physical distribution P p (h,t) 
instead of the comoving one P c (h,t). As we already mentioned, P p (h,t) describes the prob- 
ability of finding the Higgs field value h at a given time in a physical volume and takes 
into account the exponential growth of the volume of the various inflationary patches. This 
amounts to multiplying the comoving probability by exp(3iV), where N is the total number 
of e-folds. Putting aside the technical problems related to the normalizability of P p , we can 
express the physical survival probability Pa,p as the product of the exponential volume factor 
with Pa given in eq. fl36|) 



P^ : p cx exp 



+ 3N 



(41) 



16(p + 2)A 2 

From this expression it is clear that, when Hf > A8(p + 2)A 2 , the term in the exponential 
accounting for the dynamics of the Higgs fluctuations dominates over the one due to the 
physical volume expansion. Also, when ( and N are used as independent scanning variables, 
see eq. (1371) . the Higgs term exp[— a(N/N*) 1+p / 2 ] has a steeper dependence on N than the 
volume term exp(3iV), when N grows. This implies that even the physical survival prob- 
ability is pushed towards small values of the number of e-folds N and our results are not 
modified by using P p rather than P c . Because we have found that statistics prefer small iV 
(or equivalently small (pi), we can avoid problems related to the computation of the physical 
probability when this is maximized for large values of the inflaton field. Indeed, for val- 
ues of the inflaton field (p > /i(M p / fi) 6 ^ p+2 \ the inflaton evolution is quantum rather than 
classical giving rise to eternal inflation [25]. The computation of the physical probability to 
get a certain value of the cosmological perturbation ( in a physical volume would become 
therefore a complicated task as the Brownian motion of the inflaton field described by the 
Fokker-Planck equation has to be taken into account. In ref. [28] it is argued that the vol- 
ume factor exp(3iV) statistically favors small values of (. This conclusion is based on the 
requirement that the inflaton evolution remains classical, and this determines a maximum 
value of the number of e-folds N max ~ (M p //i) 2 ( 4 ~ p )/(P +2 ) ~ (- 4 /(p+ 2 ), which selects small 
(. This requirement is not well motivated on physical grounds and the quantum regime 
might modify the parametric dependence of the probability on the perturbation (. On the 
contrary, in our case the Higgs fluctuations enhance the probability of domains with small 
values of (pi, away from the region dominated by quantum evolution. 

So far, we have discussed the case in which inflation is driven by an inflaton field with a 
polynomial potential. There are two other classes of inflationary models, called small-field 
and hybrid models [2T] , where the potential is made of a constant vacuum energy plus a 
field-dependent term responsible for the slow-roll. In such a case the connection between the 
Hubble rate and the cosmological perturbation is lost. For instance, in a model with potential 
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Figure 5: Upper bound on -^/j^jif as a function of obtained by requiring that the 
electroweak vacuum is stable. Here £ is the negative coupling between the Higgs bilinear and 
the scalar curvature, and H is the maximal Hubble rate during inflation. 

V(<f>) — Vo + m 2 2 /2, the Hubble rate during inflation is roughly if* ~ (C0//A4), where 0/ 
is the value of the inflaton field at which inflation stops. Therefore, one can have a large 
probability of having domains with the correct amount of perturbation, but values of the 
Hubble rate easily smaller than the instability point A by simply choosing 0j small enough. 
However, as will be discussed in the next section, if a future measurement of gravity waves 
through the B-mode of the CMB polarization indicates that the Hubble rate is sizeable, 
even within these classes of models, one would have to conclude that we live in an unlikely 
Universe. 

Finally, we study the case in which we relax the assumption that the Higgs is minimally 
coupled to gravity. Indeed, the Higgs field might be coupled nonminimally to gravity through 
a term in the Lagrangian £Rh 2 /2, where £ is a numerical coefficient and R is the Ricci 
scalar. During inflation this term gives a contribution to the Higgs mass equal to 12C,H 2 
(since R = —12H 2 ). If < £ 10 _1 , the Higgs field is practically massless during inflation 
and our conclusions do not change. However, when |£| is sizable, our results are drastically 
modified. If £ > 10 -1 , the electroweak vacuum is stable and the Higgs fluctuations are 
damped exponentially because of the large positive effective mass. However, if £ < (as in 
the case of conformal coupling, where £ = —1/6), the electroweak vacuum is destabilized 
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by a large tachyonic mass even without inflationary fluctuations. In fig. [5] we show the 
upper bound on y/\£\H, in the case of negative £, obtained by imposing that the electroweak 
vacuum is not destabilized by the large negative mass squared induced by the £,Rh 2 /2 term. 
Here H has to be understood as the maximal Hubble rate during inflation. Note that this 
bound does not depend on the statistical considerations developed in this section, since it 
follows directly from the stability of the electroweak vacuum in our patch of the Universe. 

One could also envision the case in which there is a direct coupling between the inflaton 
and the Higgs fields, say of the form <p 2 h 2 . During inflation, this can play the role of a Higgs 
mass term, which can stabilize or destabilize the electroweak vacuum, depending on the 
sign of the corresponding coupling constant, in analogy with the case of the Higgs- curvature 
interaction. 



6 Implications for gravity wave signals 

The considerations presented in the previous sections can have direct observational conse- 
quences once we realize that the Hubble rate parametrizes the amount of tensor perturba- 
tions during inflation. During the inflationary epoch, tensor perturbations, as for any other 
massless scalar field, are quantum-mechanically generated. They can give rise to -B-modes of 
polarization of the CMB radiation through Thomson scatterings of the CMB photons off free 
electrons at last scattering [29J. The amplitude of the S-modes depends on the amplitude 
of the gravity waves generated during inflation, which in turn depends on the energy scale 
at which inflation occured. The tensor-to-scalar power ratio is often defined as T / S and is 
given by T/S ~ (if*/6.6 x 10 13 GeV) 2 . Current CMB anisotropy data impose the upper 
bound T/S < 0.6 [9j[30j. 

The possibility of detecting gravity waves from inflation via 5-modes is currently being 
considered by a number of ground, balloon and space based experiments. The decomposition 
of the CMB polarization into E- and S-modes requires a full sky data coverage and, as 
such, is limited by the foreground contaminations. The latter introduce a mixing of the 
E polarization into B with the corresponding cosmic variance limitation. Achieving levels 
below T/S = 10" 4 , or < 7 x 10 11 GeV, requires observing 70 % of the sky with dust 
emission at 0.01 % level |3ljf. 

Detection of the tensor mode would therefore imply that the value of the Hubble rate 

3 Gravitational lensing also contaminates the tensor signal by converting the dominant E polarization 
into B polarization. Cleaning this contamination by reconstructing the lensing potential from CMB itself, 
one can achieve values of T/S as small as 10 -6 [32] . 
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Figure 6: The probability Pa versus the Higgs mass normalized with its value at for 
the Hubble rate = 10 13 GeV (upper curves) or = 5 x 10 13 GeV (lower curves), 
N — iV* = 60 and for three different values of the top mass for a model of inflation with 
constant Hubble rate. 

during inflation is larger than about 10 12 GeV. This will have important implications for the 
Higgs mass within all inflationary models. In the vast class of models where the potential is 
dominated by a constant vacuum energy and the Hubble rate H* during inflation is roughly 
constant (including the small-field and the hybrid models of inflation) the regions which 
survived the Higgs fluctuations are simply distributed as 

P A oce~i^. (42) 

This probability is dominated by the smallest value of the number of e-folds and by the 
smallest value of the Hubble rate. Along the lines of the previous section, we conclude 
that a successful detection of gravity waves with a value of if* larger than about 10 12 GeV, 
in conjunction with the discovery of a light Higgs, will indicate that we live in a highly 
improbable Universe. In fig. [6] we plot the probability Pa versus the Higgs mass normalized 
with its value at = M£ for the Hubble rate ii* = 10 13 GeV (upper curves), iV = iV* = 60 
and for three different values of the top mass. We remind the reader that M£ is the upper 
metastability limit of the Higgs mass given by eq. (j2J) and shown by the upper curves in 
fig. 2. 

Similar results are also valid in the case of large-field models of inflation, with the only 
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modification that Pa is given by eq. (136]) . Therefore, in fig. [6] we also plot the probability Pa 
versus the Higgs mass normalized with its value at Mh = for the Hubble rate if* = 5 x 10 13 
GeV (lower curves), N = AT* = 60 and for three different values of the top mass. This 
value of P* corresponds to the Hubble rate necessary to reproduce the observed curvature 
perturbations for p = 2, see eq. fl35|) . 

In conclusion, our Universe becomes exponentially unlikely if measurements of the Higgs 
mass and of the tensor-to-scalar power ratio T/S violate the condition 



where the value of A has to be inferred from the measured value of the Higgs mass from 
fig. 1. The violation of eq. ( H3l) could then be interpreted as an indication for the existence of 
new physics beyond the SM, which affects the extrapolation of the Higgs potential to large 
energy scales. 

7 Conclusions 

In this paper we have investigated some possible cosmological implications of the Higgs mass 
measurement. If the LHC discovers a Higgs with mass in the metastability window shown 
in fig. |2l and does not find direct evidence for other new physics, then there is a concrete 
possibility that we live in a metastable state. If we assume that the pure SM is valid up 
to very large energy scales, then stability against field fluctuations tests properties of the 
early Universe. We have revisited the known considerations about thermal fluctuations, 
interpreting the result as an upper limit on the reheating temperature Trh after inflation. 
The bound is summarized in fig. [31 

We have also discussed the possibility that the inflationary vacuum fluctuations desta- 
bilize the electroweak vacuum. If the Hubble rate is large enough during inflation and the 
Higgs mass is light, then the danger exists that the classical value of the Higgs field is pushed 
above its instability point causing the collapse of the corresponding inflating domain. This 
does not necessarily pose a problem, since all the surviving domains will be characterized 
by the correct electroweak vacuum. However, interesting probabilistic conclusions can be 
reached under the assumption of a "landscape" scenario in which the inflationary parame- 
ters take different values in different patches of the Universe. In the context of large-field 
models of inflation we have argued that, among the surviving regions, only an exponentially 
tiny fraction of them will be characterized by the correct amount of cosmological perturba- 




(43) 
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tions. If the Higgs mass is found below (120 — 130) GeV, either we live in a very special, 
and exponentially unlikely, domain or new physics must exist below the scale A. 

Moreover, our considerations can be directly related to the amount of primordial gravity 
waves which can be measured through their imprint on the CMB. The discoveries of a light 
Higgs boson and of tensor modes would be mutually conflicting (at least in a probabilistic 
sense) if the condition in eq. ( H3|) is not satisfied. Again, this evidence could be interpreted as 
an indication of new physics modifying the extrapolation of the Higgs potential to large field 
values. This result is valid, independently of the particular inflationary model considered. 

Finally, a measurement of the Higgs mass below about 130 GeV provides a direct bound, 
shown in fig. [51 on the quantity £,H 2 , where £ is the (negative) coupling between the Higgs 
bilinear and the curvature, and H is the maximal Hubble rate during inflation. This bound 
becomes particularly interesting in case of detection of primordial gravity waves, which 
provide a measurement of H. Note that the bound in fig. [5] does not depend on any statistical 
consideration of parameter scanning. 
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Appendix 

In this appendix we present the relations between running couplings and pole masses for the 
top quark and the Higgs, and the RG equations used in our numerical calculation. We work 
at the next-to-leading order and use the MS renormalization scheme. We define 
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The corrections to the top-quark mass are given by 
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In the QCD part we include up to two-loop corrections using 
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where Nf = 5. Although the two-loop effect is beyond the precision of our computation, we 
include it because it is known to be large. The two-loop mixed electroweak-QCD [35J and 
the three-loop QCD corrections are also known [36], but give only a small effect. 



The electroweak part, including subleading corrections, is well approximated by 
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300 GeV 175 GeV 

where r = M|/(4M 2 ). The formula above is valid for r < 1, the case of interest for us. 
For r > 1 one has to replace (1/r — l) 3 / 2 arccos yjr by (1 — l/r^^arccosh-y/r. The sign of 
6.90 x 10~ 3 corrects [M] the misprint of ref. 



For the relation ([44]) between the Higgs pole mass and the quartic coupling A we use [39] 
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with f = Ml /Ml, s 
Z(z) = 



sin 2 ^, c 2 



cos 2 9w and 9w the Weinberg angle. In addition, 
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We also collect here the 2-loop renormalization-group equations that describe the evolu- 
tion of the SM couplings and Higgs wave-function renormalization at scales beyond M t [40J. 
We keep only the gauge couplings, the top Yukawa coupling and the scalar quartic coupling. 
The 2-loop RG equations for the gauge couplings ^ = {g',g,g s } are 



(54) 



where t = InQ, Q is the renormalization scale, k = l/(167r 2 ) and 

199/18 9/2 44/3 
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For the top Yukawa coupling 
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The RG equation for the Higgs quartic coupling is 
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Finally, the RG equation for the Higgs field wave function renormalization is 
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